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ABSTRACT 


Conformal mapping has been applied to determine the shape of two-dimensional 
solidified layers formed on a cold surface immersed in a flowing warm liquid. The 
frozen region is represented by a rectangle in a potential plane, and the rectangle is 
mapped into the physical plane to determine the frozen boundary configuration. The 
method is demonstrated for solidification on a cold plate of finite width that is insulated 
along its edges. 



CONFORMAL MAPPING FOR STEADY TWO-DIMENSIONAL 
SOLIDIFICATION ON A COLD SURFACE IN FLOWING LIQUID 

by Robert Siegel 
Lewis Research Center 

SUMMARY 

A method has been developed by which conformal mapping is applied to determine 
the shape of two-dimensional solidified layers formed on a cold surface immersed in 
flowing warm liquid. The liquid supplies heat by convection to the frozen interface, and 
the shape of the frozen layer adjusts so that this energy can be conducted through the 
layer to the cold surface. The interface shape is obtained by mapping the solidified 
region into a potential plane and by determining the conformal transformations between 
the potential and physical planes. The heat flow through the frozen material is also 
found. The method is applied specifically to solidification on a cold plate of finite width 
that is insulated on the edges. 


INTRODUCTION 

A method is described for computing steady-state frozen-layer configurations 
formed by solidification on a cooled plate immersed in a flowing liquid. The liquid 
is warmer than the freezing temperature, and consequently heat is transferred by con- 
vection to the frozen-layer - liquid interface. This energy is conducted through the layer 
to the cooled surface. The frozen layer will grow until the energy convected to the inter- 
face is in balance with that conducted through the layer. A steady -state frozen layer will 
thus be achieved. 

A relatively small amount of analysis has been done on freezing with either a con- 
vective boundary condition or in a two-dimensional geometry. The convective boundary 
condition was considered in reference 1, where a one-dimensional transient solution was 
obtained for solidification of a liquid flowing over a flat plate. An example of two- 
dimensional solidification is reported in reference 2, which deals with freezing inside a 
square prism with the liquid always at the fusion temperature; for this condition there 
is no heat convected to the frozen interface. The time-varying frozen-layer configur- 


ations were determined approximately. References 3 and 4 considered steady and tran- 
sient solidified layers in flow channels that are symmetric in cross section, being either 
a circular tube or a gap between parallel plates. The layers were two-dimensional since 
they varied in thickness along the channel length. However, axial conduction within the 
solidified layer was neglected so that the solidification portion of the analysis was one- 
dimensional at each axial location along the channel. The present analysis will treat only 
a steady-state situation, but it will deal with the complexities of a convective boundary 
condition and two -dimensionality. 

In the present problem the free-stream temperature of the flowing liquid and the 
convective heat-transfer coefficient at the solidified-layer - liquid interface are assumed 
constant. Then, since the freezing temperature is also constant, the convective heat 
flux supplied to the frozen interface is constant. The interface is thus one of constant 
heat flux and constant temperature; the shape of the interface will form to meet these 
constraints. 

A conformal mapping procedure is developed to determine the shape of the free 
boundary of the solidified layer. The boundary conditions on the frozen layer permit 
mapping it into a region of a potential plane. Conformal transformations are used to 
determine the frozen-layer shape from the temperature distribution represented by the 
region in the potential plane. The heat flow through the frozen layer can then be found. 
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dimensionless half width of plate, a/ y 
half width of plate 
intermediate parameter in mapping 
constants of integration 

heat transfer coefficient from flowing liquid to frozen interface 
elliptic integral 

thermal conductivity of solidified material 
length of plate 

dimensionless outward normal, n/y 
outward normal to interface 
heat flow rate through frozen layer 
heat flux 
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dimensionless frozen-layer - liquid interface, s/y 

frozen-layer - liquid interface 

dimensionless temperature, (t - t w )/(tj - t ) 

temperature 

freezing temperature 

liquid temperature 

surface temperature of cold plate 

intermediate mapping plane, £ + i-q 

intermediate mapping plane, a + i£ 

analytic function, -T + iif/ 

dimensionless coordinates x/y, y/y 

Cartesian coordinates in physical plane 

thickness of thin layer 

dimensionless physical plane, X + iY 

physical plane, x + iy 

coordinates of V plane 

length scale parameter — — — 

quantity defined as - — + i 

ax 3Y 

argument of f plane 
coordinates of u plane 
dummy integration variable 
coordinate in W plane 
mapping plane, log | £ | + id 


ANALYSIS 

This analysis will develop an analytical method which applies conformal mapping to 
determine the shape of two-dimensional steady -state frozen regions formed in a flowing 
liquid. The method will be applied to the specific geometry of freezing on a constant - 
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Flowing liquid, tj,h 



Figure 1. -Two-dimensional frozen layer formed on cold plate 
immersed in flowing liquid. 


temperature plate of finite width as shown in figure 1. The plate, of width 2a, is insulated 
on its two sides; the length of the plate L is large, so that the configuration can be con- 
sidered two dimensional. A liquid with a bulk temperature t^ that is higher than the 
freezing temperature tj is flowing in the direction along the plate length. The flow pro- 
vides a constant heat-transfer coefficient h between the liquid and the surface of the 
frozen layer at tj. 

Because h, , and tj are all constant, the heat flux transferred to the surface of 
the frozen layer is uniform over that surface. Thus, the boundary conditions on the 
frozen -layer - liquid interface are both uniform heat flux and uniform temperature: 

' = h(t ? - t f ) (la) 

s 

t(s) = t f (lb) 

The shape of the frozen-layer - liquid interface will adjust until the heat transferred to it 
by the flowing liquid is exactly balanced by the heat conducted through the frozen layer to 
the cooled plate. 

When utilizing conformal mapping it is best to define dimensionless variables so 
that the dimensionless temperature and its normal derivative are unity at the frozen 
interface. Thus, the variable is defined as 


q(s) = k — 

3n 


T = 


t - t. 


w 
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Then, at the plate surface 


T = 0 (-a < x < a; y = 0) 

and at the frozen-layer - liquid interface 


(2a) 


T(s) = 1 


The boundary -condition equation (la) becomes 


0T 

3n 


b h ■ *f 


s k ‘l - ‘w 


A length scale parameter y is defined as 


i t. - t 
k f w 

y = 

h t, -t. 


(2b) 


so that the temperature derivative can also be normalized to unity: 



(2c) 


With all lengths nondimensioned by y (i. e. , X = x/y, Y = y/y, A = a/y), the frozen 
layer with its associated boundary conditions is shown in figure 2 in the dimensionless 
physical plane Z = X + iY. 


iY 



Figure 2. - Frozen layer in dimensionless physical plane, Z - X + iY. 
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The determination of the shape of the interface between the frozen layer and liquid 
is similar to the computation of the boundary of a free inviscid jet flowing into a constant- 
pressure environment. The velocity along the jet boundary is a constant as specified by 
Bernoulli’s equation, but the jet boundary configuration is unknown. In the present case 
the frozen-layer - liquid interface geometry is again unknown, and the boundary conditions 
along the frozen-layer - liquid interface are known constants. 

The conformal mapping procedure used for the freezing problem is based on the fact 
that the temperature distribution within the frozen layer is a solution to the Laplace heat 
conduction equation. Since the frozen region is simply connected, the temperature dis- 
tribution can be taken as the real part of an analytic function W. In the case of heat con- 
duction, the positive flow of heat is in the direction of the negative temperature gradient. 
As a result the potential function (real part of W) is -T, and the analytic function can be 
written as 


W = -T + ty 


( 3 ) 


The lines of constant T and constant xf/ form an orthogonal net. By using conformal 
transformations between the W plane and the physical plane, the transformed functions 
will always be analytic and, hence, in the physical plane the temperatures will satisfy 
the Laplace equation. Since a function which is analytic in a given region is completely 
determined by its boundary values, it is only necessary to require that this mapping 
satisfy the boundary conditions. 

The isotherms and adiabatic lines bounding the frozen layer are shown in figure 3, 
where the layer becomes a rectangle. Vertical lines in the W plane are isotherms, and 
horizontal lines are heat flow lines that are normal to the isotherms. Hence, for the 
present problem the segments 2-3 and 5-6 in the Z plane become horizontal lines in the 
W plane since along these lines 3T/ 3Y = 0 and they are consequently normal to the 
isotherms. The height of the rectangle is not known but will be determined later to 
obtain the heat flow through the frozen layer. 



--T 


Figure 3. - Potential plane, W = -T + i 



If the mapping between the W and Z planes is known, then the line 2-6 in the W 
plane can be transformed back to the physical plane to determine the interface configur- 
ation. After the mapping is completed, the heat flow through the frozen layer can be 
computed. 

To find the relation between the W and Z planes, the condition is utilized that the 
value of the derivative of an analytic function is independent of direction. Consequently, 

dW = _9T + il£ (4 

dz ax ax 

Because the temperature distribution within the steady-state frozen layer is governed by 
the Laplace equation, the Cauchy-Riemann equations can be used to further relate deri- 
vatives of x{/ and T. Using the relation 


d\j/ _ 3T 
ax 3Y 


equation (4) can be written as 


dW _ 3T . 3T 
dZ 3X 3Y 


Now a new quantity £ is defined as 


Then equation (6) becomes 


£ = 


3T i 3T 
ax 3Y 



After integrating, 


( 7 ) 


Z = J - dW + C x (8) 

where C ^ is a constant of integration. Equation (8) shows that the connection between 
the physical plane and the potential plane depends on the function £ which must be related 
to W before the integration can be performed. 
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Figure4. - Temperature derivative plane, £=-|^- + i|y. 



To examine the relation between £ and W, a temperature derivative plane is con- 
structed, yielding figure 4 for the present problem. The interface 2-1-7-6 has 
3T/ 3N = 1 all along it since the interface is at constant temperature and the heat flow is 
normal to the constant temperature lines. As a result the interface transforms to a unit 
semicircle on the £ plane. Along 2-3 and 5-6, 3T/ 3Y = 0 because the boundary is 
perfectly insulated, but 3T/ 3X varies from +1 and -1 at points 2 and 6 to +<*> and -«> at 
the corners 3 and 5 where the heat flow lines have infinite curvature. Along 3-4-5, 

3T/ 3X = 0 since the boundary is at constant temperature. The 3T/ 3Y varies from <*> 
at the corners 3 and 5 to a value at 4 that is unknown at this point in the analysis. 

To map the region in the £ plane into the W plane and thereby determine the 
relation between £ and W, a logarithmic transformation is first applied: 

w = log £ = log | £ | + i# (9) 

where 6 is the argument of £. This maps the semicircle into a straight line and yields 
the configuration in figure 5. 


ie 
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Figure 5. - Intermediate u plane, u ■ log |£| +i6. 
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The region in figure 5 is in the form of a generalized rectangle which can be mapped 
into the positive half plane by using a Schwarz -Christoff el transformation. This trans- 
formation has sufficient degrees of freedom so that the points 3 and 5 can be mapped to -1 
and +1 on the real axis. Points 2 and 6 are at -b and +b, where |b | < 1 and is a quantity 
that is to be determined. The transformation is 


dio 

du 



( 10 ) 


and the configuration in the u plane is shown in figure 6. 


it) 



Figure 6. - Intermediate u plane, u = £ + ly. 


Equation (10) is integrated to yield 



To have point 1 at co = in/ 2 correspond to u = 0, the constant Cg = i7r/2. Also, as u 
approaches b and -b, o> approaches 0 and in, respectively; these conditions require 



The upper half of the u plane of figure 6 is then mapped into a rectangle in the V 
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Figu re 7. - Intermediate V plane, V = a + ip. 


plane as shown in figure 7. The mapping transformation is found from the Schwarz - 
Christoffel theorem as 


dV 

du 


The final mapping from the V 
tion given by 


Vu - b Vu + b Vu - 1 Vu + 1 

plane to the W plane of figure 3 is a minus 90° 


( 12 ) 

rota- 


W = -iV 


(13) 


With the mapping relations established between f and W, equation (8) can be used to 
transform to the dimensionless physical plane by writing the integral in the form 


Z = 



1 dW dV 
£(u) dV du 


du + C j 


(14) 


To find £(u), equation (9) is used in the form f = e w . Substituting ui from equation (11) 
and noting that = i gives 



( 15 ) 


Then with dV/du from equation (12) and dW/dV = -i, equation (14) becomes 
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Multiplying both numerator and denominator of the term in parentheses by 



This can be further simplified to 



Equation (17) is integrated to obtain 



where there is a branch cut between -b and +b along the real axis of the u plane. 

The constants are evaluated by having Z = 0 at u = ±«o and Z = +A and -A at 
u = -1 and +1, respectively. This gives 



Then, 
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(19) 


£ = 1 = i (log Vu 2 -b 2 - Vi-b 2 . 

A 2 log Vl -b 2 \ Vu 2 - b 2 + Vl - b 2 U + > 

The transformed frozen-layer - liquid interface extends between | = +b and £ = -b 
on the real axis of the u plane in figure 6. If £ is introduced for u in equation (19), 
where 1 1 [ < b, then the real and imaginary parts can be computed to obtain the x and 
y coordinates of the frozen-layer - liquid interface: 


x 

a 


s 


2 log 


■ - log 

Vl - b 2 



y 
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s 


-1 


tan 


-1 


log 


VTb 2 



(20a) 


(20b) 


where -b < £ < b. 

The quantity b must now be related to the imposed physical conditions. It will be 
found that b is a function of only the single parameter [(t f - t w )/(t, -t f )](k/ha). This will 
also lead to a relation between this parameter and the total amount of heat being trans- 
ferred through the frozen layer. 

The heat flow through the frozen layer and into the cold plate can be written in terms 
of the temperature gradient at the plate: 


Q 



( 21 ) 


where the integration limits are the numbered positions in figure 2. Equation (21) is 
placed in dimensionless form: 


Q 

2kL(t f - t w ) 



(22) 


By using equation (5), equation (22) can be expressed in terms of the stream function 
and then integrated: 
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Q 

2kL(t f - t w ) 



MdX = tf/(3) - i//(4) 
8X 


From figure 3, rf/( 3) = ^(2) and \f/( 4) = i//(l), so that 


Q 

2kL(t f - t w ) 


= i//(2) - ^(1) 


( 23 ) 


(24) 


For the problem being considered here equation (24) is more convenient than equation (23) 
because points 2 and 1 are at finite values of u in figure 6, while point 4 is at 
u = ±oo + iO. This will make the integration limits finite in equation (25). 

To find \p(2) - i//(l) and thereby evaluate Q, equation (3) is used: 


■M2) - f(l) 



ildu 
dV du 


(25) 


where the integration is performed for u along the real axis in figure 6. Then, using 
equations (12) and (13) and noting that u = £ on the real axis, 


Q 

2kL(t f - t w ) 




(26) 


The constant was previously evaluated as 

Ab 


c 4 = 


ha h " h 


I o k t- - t. 

log Vl -b 2 f w log Ml -b 


f L w ^2 


Inserting equation (27) into equation (26) yields 


Q 


2kL(t f - 


r - b 

ha h ' b 1 dj_ 

k h - v log Vi _ b 2 J o Vb 2 - ? VT7 


(27) 


(28) 
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A second relation for can be found by noting that in figure 3 the size of the 
rectangle must be such that T(2) - T(3) = 1. Then in a manner similar to the derivation 
of equations (25) and (26) 



Solving for and inserting into equation (26) yields 



(29) 


By substituting equation (29) for Q into equation (28) the result can be rearranged into 


k_ *f ~ *w 
ha tj - tj 


log 






V? 2 - b 2 Vi - 7 


(30) 


The integrals in equations (29) and (30) are elliptic integrals that can be placed in 
the following forms by using trigonometric substitutions (ref. 5): 


Q 


K(b) 


2kL(t f - *-> k(VT7) 


(31) 


k ~ *w 
ha tj - tj 


K 


log 




(VTT^) 


(32) 
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where 



Equations (20), (31), and (32) constitute the solution for the frozen layer being con- 
sidered here. The physical quantities are related implicitly through the quantity b. 
When a value is chosen for b between 0 and 1, equations (20) give the configuration of 
the frozen layer, equation (31) gives the heat flow through the layer, and equation (32) 
gives the magnitude of the physical parameter (k/ha)|(tj - ^ w )/(t^ - t^)J associated 
with the layer. 


RESULTS AND DISCUSSION 


By using equation (32) the values of b corresponding to various 
(k/ha)[(t f - t w )/(tj - t f j] can be found. Then these b values are used in equations (20) 
to obtain the coordinates of the frozen-layer - liquid interface. The frozen-layer con- 
figuration is shown in figure 8. 

Thin layers result when the parameter (k/ha)[(tj - t w )/ (t^ - t^)j is small (fig. 8(a)). 
This would be expected because a small value of the parameter occurs when the cooling 
of the plate is small and there is a large convective energy flux supplied to the frozen 
interface. The frozen layer forms to the proper configuration so that the convective flux 
can be conducted through it to the cold plate. For very thin layers, the layer approaches 
a constant thickness, and the heat flux through the layer is one dimensional. Then, let- 
ting the one-dimensional thickness be y^ n , 


h(t z 



k(t f - t w ) 

^thin 
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Solving for the dimensionless thickness, 


y thin _ k *f ~ *w 
a ha t, - tj 


The thickness in figure 8(a) thus approaches the value of the parameter 
(k/ha)[(t f - t w )/(t z - t f )] as the parameter becomes small. 

For large values of (k/ha) j(tj - t w >/% - 1 £ )J, as shown in figure 8(b), the frozen 
layer becomes quite thick and approaches a circular shape. In the limit of very large 

the width of the plate is small relative to the frozen-layer 


(k/ha)[(t 


y/(t; 


l 


t >] 


thickness, and the configuration approaches a two-dimensional sink with the heat flowing 


radially through the layer to the cold plate. 

By inserting various b values into equations (31) and (32) the dimensionless heat 
flow through the layer Q/2kL(tj - t w ) can be evaluated in terms of the parameter 
(k/ha)|(t£ - t w )/(t^ - tj)J. This is shown in figure 9. As expected, the heat flow decreases 



continuously as the frozen-layer thickness increases. For very thin layers the heat flow 
becomes one -dimensional and flows parallel to the y-axis. The surface area of the fro- 
zen interface approaches 2aL. Then the heat flow is Q = h(t^ - t^)2aL, or, in dimen- 
sionless form, 
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Q _ ha H ~ *f 
2kL(t f - t w ) " k t f - t w 


CONCLUDING REMARKS 

Conformal mapping has been applied for determining two-dimensional steady-state 
configurations of a solidified layer formed on a cooled plate immersed in a flowing liquid 
The frozen layer is mapped into a potential plane W = -T + ixjy where the negative of the 
dimensionless temperature is the potential function for heat flow. The potential plane is 
related to the physical plane by means of a function involving temperature derivatives: 


£ 


9T i 0T 
aX 3Y 


The conditions at the frbzen-layer boundary are used to determine the form of the £ 
function. Then conformal mapping is applied to relate £ to W, and the physical config- 
uration is found from the integral 

Z = f -J— dW 

J £(W) 

The method was demonstrated specifically in the analysis of solidification on a 
cooled plate of finite width. An example for further application would be to determine the 
shape of solidified material inside noncircular liquid flow channels with cold walls. These 
types of channels are found in heat exchangers that are cooled cryogenically. 


Lewis Research Center, 
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